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Auxiliary- 
Library Routine riA 1 - 257 



Floating Point Simpson's Rule Integration (SADOI Only 

Closed 

J+l 

The truncation error in Simpson's Rule in taking 2 

steps each of length h is: 

ilih 

90 
where i is some point in the interval. 
[(52s+ a) + 75 + a] m.s. where s is the number of 
steps and A is the time in m.s. required by the 
auxiliary. 
An approximate value for the integral 

I =J a f(x) dx 
is obtained by 

I =^{f(xQ) + f(xj + 4[f(x^) + f(x^) + .., 

+ f (Xjj.-l)] + 2[f (x^) + f (x^^) + ... 



There must be an even number of 



+ f 

where h = x. - x. ^ . 

iirteirvala, n, givea hy an odd number, n + 1, of points 

Xq, x^, ..., x^. The minimal number of points is 3, for 

which I = 5 [ f(xo) + ^f(x^) + fCx^)]. 

The size of the interval h can be changed at any 
point with n even; thus, for example we may have for 
two intervals 

I = ^ [f(x ) + i^f(xj + 2f(xJ + kf(x.) 
-^ h^ ^ 

+ f (x^)] + j-[f (xj^) + ^f (x^) + 2f (x^) 

+ iff(x^) + ... + f(x^)]. 
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EMRY: 



EXIT: 



AUXILIARY: 



Entry is made from fixed point in the standard vay with 
control transferred to symbolic address (SIMP), i.e. 

50 p 



p + 1 



26 (SIMP) 



Upon entry the following parameters must be stored in 
the following order at symbolic location (l): 

""l^ ^1' ^01* '^' \' ^02^ ••• ^ V V ""ok* ° 
with the zero following the last group to denote the 
termination, where n (a positive, even, fixed point 
integer, ie. n^ x 2-'9) = number of steps to be taken 
with this step length, h^ (a floating point number) 
= length of step, x^^ ( a floating point number) 
- starting point of integration with this length. 

Hence the step length can be changed without leaving the 
routine . 

Control is returned, in fixed point, to the right hand 
order in p + 1 with the value of the integral in the 
floating accumulator. 

The auxiliary subroutine for computing (or otherwise 
determining), f (x) must be located at symbolic location 
(AUX). It is entered from EA 1 in fixed point using 
an 8j order, with x in the floating accumulator. It 
should begin as follows: 



(AUX) 



L5 2S4 
^0 (SAVE) 
52 l(AUX) 
50 l(AUX) 
2E~^ 



Waste 

Enter Al with x in the 
floating accumulator 



The auxiliary should place the value of f(x) in the ■ 
floating accumulator and exit should be made as follows! 

L5 (SAVE) 

40 2S4 



26 29S4 
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ASSEMBLY : 



The floating ac::umulator must be in the locations 
specified by 35 and IS^. Lihrary routine A 1 must 
begin at the location specified by S 4. If the user 
desires other of the floating point subroutines he 
will find it most expedient to use routine A 6 by 
calling it from the drum. 

Routine EA 1 uses b -boxes and 1 of A 1 and the 
following symbolic addresses: 

(1), (J), (K), (N), (S) and (SIMP). 
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Location 



(SIMP) 



5 

k 

5 
6 

T 
8 

9 

10 
11 
12 

13 
li^ 

15 
16 



ORDER 

00 K 
K5 F 
k2 5L 
i^l(S) 
L5 12L 

Li^ 53L 
if2 7L 

k2 8l 

LI (L) 

40 F 

32 F 

LI F 

10 IF 

k2 15L 

L5 1(L) 
i^0 (K) 

L5 2(L) 

ko (n) 

^1 (J) 
^1 1(J) 
22 lifL 
L5 ^L 
Ih 52L 
26 2L 
00 (l) 
00 F 
00 F 
00 F 
50 14l 
26 Sk 
IK [^i/2]F 
85 (N) 
8j (AUX) 



NOTES 
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by 2L 
link 
^i/2 
by 5L 
by kL 



27 L 



J 



by 7L 




Clear grand sum 



A = - n. 

+ -> done 
A = n 

A = h. 



A = X, 
Set X 



Oi 




Enter F.P. 



Set loop 
Get f (xq) 



LOCATION 


ORDER 


NOTES PAGE 2 


EA 1 


IT 


8S (1) 








18 


OK 2F 
85 (N) 
^ (K) 


.51 L 


Set binary switch at 20L 
X -> X + h 




. 29 L 








19 


8S (N) 
8j (AUX) 








20 


05 28l 
15 50L 




Transfer every other time 
Transfer 1/2 - 1 times 




21 


^ (1) 
8S (1) 








22 


8k Uf 










87 (J) 








25 


^ (1) 
8S (1) 








2k 


8K 2F 
87 l(J) 








25 


^ (1) 
86 5^L 




t3 




26 


87 (K) 
m (S) 




Multiply by h^ 
Add to grand siim 




27 


8S (S) 
8J IIL 








28 


^ (J) 

8s (J) 


yfrom 20L 










29 


8k F 
85 i8l 




Set F = 




50 


^ 1(J) 
8s 1(J) 


, from 20L 




1 






51 


8k f 
82 17L 




Set F = 




52 


00 F 

00 5F 








1 55 


00 F 








1 
j 


00 IF 








5^ 


26 1658F 




= 5 in F.P, 




J - ^..^ 


66 172 9F 






. 



56 

57 

59 
40 




(1) 00 F 

00 F 
(j) 00 F 

00 F 
00 F 
00 F 

(K) 00 F 
00 F 

(k) 00 F 
00 F 

(S) 00 F 
00 F 
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